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摘 要 :通过 溃 坝 案例 的 现场 调查 和 水 楼 模型 试验 ,建立 了 面板 砂砾 石 坝 漫 顶 溃 坝 数学 模型 , 以 模 
拟 水 流 特性 、 砂 砾石 料 以 及 混凝土 面板 三 者 之 间 的 相互 作用 。 该 模型 主要 由 3 部 分 组 成 :人 考虑 湾 
涡 水 流 和 砂砾 石料 之 间 的 冲 蚀 特 性 ,引入 RNG kk-e 潢 流 模 型 模拟 洲 涡 水 流 的 反 向 侵蚀 ;@) 通 过 泥水 
输 移 公式 和 流体 体积 法 (VOF 法 ) 追 踪 砂 粒 与 水 流 侵蚀 交界 面 , 并 考虑 了 砂砾 石料 的 孔 阶 特性 ;@) 
基于 弯 适 平衡 法 判定 混凝土 面板 在 自重 和 水 荷载 作用 下 的 破坏 过 程 。 通 过 水 楷模 型 试验 对 数学 模 
型 进行 验证 ,结果 表明 ,建立 的 数学 模型 能 够 准确 地 模拟 溃 决 过 程 中 湾 涡 水 流 对 大 坝 的 反 向 侵蚀 和 
面板 溃 口 的 详细 演化 过 程 ;该 模型 计算 的 溃 口 流量 、 扒 石 体 溃 口 的 发 展 、 溃 坝 历时 以 及 面板 的 折断 
长 度 ,与 实测 数据 的 相对 误差 均 小 于 15% ,验证 了 模型 的 合理 性 。 与 参数 化 模型 相 比 ,本 研究 提出 
的 演 坝 模拟 数学 模型 具有 更 详细 的 结 
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dam breaching due to overtopping 
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Abstract :Based on the field investigation of dam break cases and flume model test ,the numerical model of 
overtopping break of concrete-face sand-gravel dam ( CFSGD ) is established to simulate the interaction a- 
mong water flow characteristics , sand-gravel materials ,and concrete face slab. The model is mainly com- 
posed of three parts: (Q) Considering the erosion characteristics between the vortex water flow and sand- 
gravel materials ,the RNG k-e turbulence model is introduced to simulate the reverse erosion of the vortex 
water flow. CC) Through the sediment transport formula and fluid volume method (VOF) traces the inter- 
face between sand and water erosion , and considers the pore characteristics of sand-gravel materials. (3) 
Based on the moment equilibrium method ,the failure process of the concrete face slab under its self-weight 
and water loads is calculated. The numerical model is verified by the flume model test. The results show 
that the numerical model established in this paper can accurately simulate the reverse erosion of the vortex 


water fow on the dam during the failure process and the detailed evolution process of the face slab breach. 
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The relative error between the measured data and the calculated peak breach flow ,the development of the 


breach of the rockfill body ,the duration of the break and the length of the face slabs are all less than 


15% ,which verifies the rationality of the numerical model. Compared with the parametric model ,the nu- 


merical model of dam breach proposed in this paper provides more detailed results. 
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water flow ;face slab fracture 


面板 砂砾 石 坝 因 其 工期 短 .就 地 取材 、 对 地 形 地 
质 条 件 适应 性 强 和 大 坝 蓄 水 期 或 运行 期 的 坝 体 沉降 
变形 小 等 优点 , 逐渐 成 为 一 种 具有 竞争 力 的 坝 
型 5 。 目 前 ,国内 外 越 来 越 多 高 面板 坝 采用 砂砾 石 
料 填 筑 坝 体 。 例 如 , 坝 高 187 m 的 Aguamilpa 坝 和 坝 
高 123.5 m 的 青海 黑 泉 水库。 然而 砂砾 石料 的 抗 冲 
促 性 能 较 差 ,一 旦 因 超 标 洪水 或 库 区 滑坡 等 因素 
库 水 漫 顶 ,面板 砂砾 石 坝 将 极 有 可 能 发 生 省 坝 。 
1993 年 , 坝 高 71 m 的 沟 后 水 库 发 生 溃 坝 事故 ,造成 
房屋 良田 被 毁 和 320 人 死亡 5 。2005 年 , 坝 高 
2734 m 的 Taum Sauk 上 水 库 因 水 位 计 失 灵 造 成 库 水 
漫 斋 , 仅 25 min 内 5.3 x 105 m? 的 水 汇 入 Black 
河 回 。 为 了 减轻 或 避免 省 坝 造成 的 损失 ,开展 面板 
不 故 石 坝 漫 顶 溃 坝 数值 模拟 研究 具有 重要 的 意义 。 
EN 国内 外 学 者 针对 均 质 坝 漫 顶 溃 决 机 理 展开 深入 
5 ,得 出 坝 体 表 层 冲刷 “ 陡 坎 式 " 后 退 冲 蚀 以 
贵 口 边 坡 失 稳 是 均 质 坝 漫 顶 溃 坝 的 基本 过 程 。 然 
= 由 于 混凝土 面板 的 支撑 作用 ,面板 坝 的 漫 顶 溃 决 
析 理 与 均 质 坝 明显 不 同 。 沟 后 水 库 和 Taum Sauk 上 
水 座 溃 坝 后 的 调查 结果 表明 , 坝 体 游 涡 冲 刷 ` 面 板 往 
复 牟 断 以 及 溃 口 边 坡 失 稳 是 面板 坝 溃 坝 的 3 个 主要 
特 钨 ,如 图 1 所 示 。 基 于 沟 后 面板 砂砾 石 坝 溃 决 案 
例 , 刘 杰 等 5 通过 水 槽 模型 试验 查 清 了 沟 后 面板 砂 
砾石 坝 失 事 的 原因 ,但 未 揭示 库 水 漫 顶 后 大 坝 的 溃 
决 机 理 ; 李 雷 等 "1 采用 SINGH 提出 的 冲 蚀 公式 建 
立 了 简化 演 坝 数学 模型 ,但 未 考虑 演 口 边 坡 的 破坏 ; 
王 廷 等 "考虑 下 油水 流 中 泥 沙 含量 变化 对 溃 决 过 
程 的 影响 ,建立 了 动 床 粗 合 分 析 模型 ,但 对 面板 折断 
判定 分 析 过 于 简化 。 根 据 实际 溃 坝 案例 调查 和 离心 
模型 试验 , 陈 生 水 等 '“ .ZHONG 等 '“ "1 详细 研究 了 
砂砾 石料 的 输 移 和 混凝土 面板 的 破坏 ,提出 了 面板 
砂砾 石 坝 漫 顶 溃 坝 数学 模型 ,但 忽略 了 洲 涡 水 流 对 
坝 体 的 冲 蚀 与 淘 刷 。 
Flow-3D 软件 采用 独特 的 TruVOF 技术 和 经 典 
的 泥 沙 冲 蚀 模型 ,可 以 准确 地 计算 砂粒 的 沉积 、 冲 刨 
以 及 运动 ,并 能 够 很 好 地 模拟 溃 坝 过 程 中 复杂 的 水 


砂 流动 状态 。 本 研究 通过 对 Flow-3D 软件 进行 二 次 
开发 ,采用 弯 矩 平衡 法 对 混凝土 面板 的 折断 进行 判 
定 ,建立 一 种 能 够 反映 水 流 特 性 、 砂 砾石 料 和 混凝土 
面板 三 者 之 间 相 互 作用 的 漫 顶 演 坝 数学 模型 。 重 点 
考虑 洲 涡 水 流 冲 蚀 特性 水 砂 侵蚀 交界 面 和 混凝土 
面板 结构 破坏 特征 ,计算 结果 对 面板 砂砾 石 坝 建立 
风险 等 级 和 风险 预警 机 制 具有 重要 价值 。 


防 浪 墙 基础 裸露 


图 1 面板 砂砾 石 坝 演 坝 案例 


Fig.1 Dam breach case of concrete-face sand-gravel dam 


1 溃 坝 模拟 数学 模型 


基于 面板 砂砾 石 坝 的 溃 决 特点 ,模拟 面板 砂砾 
石 坝 漫 顶 溃 坝 过 程 应 着 重 考虑 3 个 方面 :@ 针 对流 
涡 水 流 对 坝 体 的 冲 创 与 淘 刷 ,采用 合适 的 模型 模拟 
水 流 运动 ;@ 准 确 地 模拟 大 坝 破坏 过 程 中 水 流 与 砂 
粒 的 侵蚀 交界 面 ;@ 分 析 面 板结 构 受 力 ,判定 面板 的 
折断 长 度 与 时 间 。 
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Meyer-Peter & Muller 方程 为 


1.1 RNG k-e 党 流 模 型 


面板 人 砂砾 石 坝 演 决 过 程 中 , 演 坝 水 流 在 面板 下 
部 区 域 发 生 剧烈 的 淘 刷 , 挟 砂 水 流 均 是 消 流 ,因此 采 
用 消 流 模型 模拟 泪 坝 水 流 。RNGC k-e 漠 流 模型 考虑 
了 流动 中 水 流 的 旋转 情况 ,能够 更 好 地 对 流 场 中 流 
线 弯曲 程度 进行 较 大 的 改进 ,并 且 该 模型 的 系数 是 
通过 理论 分 析 推 导出 来 ,计算 精度 高 于 其 他 经 验 模 
型 ,在 模拟 复杂 水 流 冲 刷 方面 得 到 了 广泛 的 应 用 。 
其 基本 控制 方程 ”为 


亲 动 能 方程 
a(pk) (pkui) 9 ,Ok 

a ll oe WS 

(1) 
巡 动 能 耗 散 率 。 方程 
©O ao(pe) (peui) 9 Mr \ oe 
Bt 1 ox, | lp | 
£ £ 

己 44-Cp 第 (2) 


其 用 :: 为 时 间 ; 为 亲 流 夭 性 系数 ; p= CPP/e ， 
(全 9 085 ;6 为 训 动 能 4 的 产生 项 。 


1QY 泥 沙 输 移 与 侵蚀 
©O 


CN 溃 决 初期 ,在 下 泄 水 流 作用 下 ,部 分 细小 颗粒 开 

始 性 生 运动 , 随 着 下 油水 流 强度 的 逐渐 增强 , 越 来 越 
多 的 砂粒 由 沉积 状态 转变 为 悬浮 状态 ; 当 溃 口 流量 
逐 嘻 减 小 时 , 溃 坝 水 流 的 挟 砂 能 力 降低 ,悬浮 砂粒 在 
自重 作用 下 开始 下 沉 。 整 个 溃 决 过 程 中 ,砂粒 的 运 
动 颈 态 在 悬浮 与 沉积 之 间 相 互 转化 。 而 下 泄 水 流 中 
挟 三 含量 的 变化 会 引起 水 流 密 度 与 条 性 系数 的 变 
化 ,进而 影响 侵蚀 过 程 。 为 了 使 计算 结果 更 加 准确 ， 
本 研究 定义 了 砂粒 与 水 流 的 混合 流体 ,该 混合 流体 
的 密度 为 


p= apv + (1 -ar)p。 (3) 
其 中 :a 为 砂粒 体积 分 数 ;p, 为 砂砾 石料 密度 ;p, 为 
水 流 密度 。 

由 于 无 法 计算 单个 颗粒 在 流体 中 的 运动 ,因此 
计算 砂粒 输 移 需 采用 经 验 公 式 。 本 研究 基于 泥 沙 携 
带 提升 速度 方程 "和 Meyer-Peter & Muller 方程 :201 
计算 砂粒 的 悬 移 质 运 动 与 推移 质 运动 。 

泥 沙 携带 提升 速度 方程 为 


wns = mde (0 ~ 9%) [el (Ee)d] 


(4) 


qs =B (0 00)7 | el real (5) 
其 中 ;n, 为 斜坡 垂直 方向 ; d, 为 无 量 纲 粒 径 参 数 ; 0， 
为 床 面 希 尔 兹 数 ; 0。; 为 临界 希 尔 兹 数 ;d 为 颗粒 平 
均 粒 径 ; 常 规 参 数 '””1 取 值 a, = 0.0018,B =8。 

模型 假定 江口 形状 为 梯形 ,将 演 口 的 发 展 分 为 
两 个 阶段 ,第 一 阶段 假定 省 口 在 深度 和 宽度 上 以 同 
样 速度 均匀 发 展 , 直 至 溃 口 发 展 至 底部 ;第 二 阶段 泪 
口 向 两 侧 扩展 ,直至 省 决 过 程 结 束 。 为 了 能 够 准确 
地 描述 坝 料 与 水 流 侵蚀 交界 面 ,VOF 法 2 (流体 体 
积 法 ) 是 一 种 有 效 的 方法 ,该 方法 通过 单元 内 流体 体 
积 与 单元 体积 之 比 来 追踪 自由 液 面 。 本 研究 在 VOF 
法 的 基础 上 考虑 了 砂砾 石料 的 孔隙 特性 ,单位 体积 
的 砂砾 石料 中 土 颗粒 仅 占 一 定 份额 ,其 余 为 空气 和 
水 。 因 此 ,单元 内 砂砾 石料 的 最 大 填充 应 根据 土 样 
的 孔隙 度 进行 确定 。 


1.3 面板 折断 


在 混凝土 面板 砂砾 石 坝 溃 决 过 程 中 , 堆 石 体 的 
沉 口 发 展 逐 渐变 深 变 冤 ,混凝土 面板 的 其 空 长 度 和 
宽度 也 随 着 堆 石 体 江口 的 发 展 而 增 大 , 当 混 凝 土 面 
板 无 法 承受 自重 和 水 荷载 作用 时 ,面板 发 生 断 裂 ,如 
图 2 所 示 。 一 块 混 凝 土 面板 折断 后 , 堆 石 体 溃 口 继 
续 横 向 扩展 和 下 切 ,随后 两 侧 的 混凝土 面板 相继 发 
生 折 断 。 此 外 , 由 于 省 口 的 下 切 ,每 块 混凝土 面板 可 
能 会 出 现 多 次 折断 。 本 研究 将 面板 假定 为 悬 岁 板 ， 
采用 弯 矩 平衡 法 对 混凝土 面板 的 折断 进行 判定 。 


区 


图 2 混凝土 面板 受 力 分 析 示 意图 


Fig.2 Schematic diagram of the force analysis 


SN 
2 


of concrete face slab 
面板 自重 产生 的 弯 和 矩 MM 可 表示 为 
RE pngmIOwL (6) 
2 V1 +m?i 
其 中 ;p, 为 混凝土 面板 密度 ; mi 为 上 游 坡 比 ; 6 为 面 
板 的 等 效 厚度 ; 上 为 面板 折断 的 极限 长 度 ; zi 为 面 


板 顶 部 高 程 ; z, 为 库 水 位 高 程 。 
水 荷载 作用 在 面板 上 的 弯 矩 Ms 可 表示 为 
vgw(z 2) pg 
2 Ir 
_Pw8w [Zu (zr -z)] 
/I+rm 
因此 ,作用 于 混凝土 面板 的 总 弯 算 为 
M = M,+M, (8) 
混凝土 面板 的 极限 弯 矩 Mj 可 根据 《水 工 混凝土 
结构 设计 规范 》 进 行 计算 ,面板 折断 的 判定 条 件 为 
M>M, (9) 
当 判 定 面板 发 生 折断 时 ,需要 对 折断 部 分 的 壁 
面 边界 条 件 进行 重新 设 定 ,使 其 不 再 具有 挡 水 作用 ， 
以 王 模 拟 混凝土 面板 的 折断 。 值得 注意 的 是 ,尽管 
模型 面板 与 实际 钢筋 混凝土 面板 存在 差异 ,但 是 数 
电 。 中 与 面板 性 能 相关 的 是 极限 弯 矩 Wi ,因此 在 
证 算 参 数 选取 时 , 按 模型 面板 或 钢筋 混凝土 面板 的 
实际 性 能 确定 即 可 。 
水 流 矫正 机 
摄像 类 ”测量 标杆 


M, = 人 


Z。>> Zr 


(7) 
M, 


Zs <Zr 


进 水 口 
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(a) 试验 装置 示意 图 


2 面板 砂砾 石 坝 演 坝 模型 试验 


2.1 试验 装置 


物理 模型 试验 在 图 3 所 示 的 钢 制 水 槽 中 进行 ， 
水 槽 前 段 长 3 m, 水 平 放 置 ,后 段 长 5 m, 具 有 5° 的 坡 
度 。 在 水 覃 的 折 坡 处 设置 有 水 流 矫正 栅 , 可 以 有 效 
地 削弱 上 游 水 流 对 模型 坝 的 影响 ,未 端 设 置 有 沉 砂 
池 ,以 便 模 型 料 的 回收 利用 。 本 试验 中 上 游 水 库 恒 
定 供 流量 为 5 L/s。 坝 体 上 游 坡 脚 处 布设 有 水 位 计 ， 
以 实时 监测 坝 前 水 位 动态 变化 过 程 。 坝 体 项 部 、 下 
游 1/2 坝 坡 处 以 及 坝 体 中 间 纵 剖面 上 布置 有 测量 标 
杆 ( 标 注 有 刻度 的 细 铁 杆 ,底部 连接 有 小 铁 片 , 填 和 贷 
料 被 冲 蚀 后 依靠 自重 下 落 ) ,测量 溃 口 的 扩展 过 程 以 
及 坝 体 纵 剖 面 的 演化 过 程 ,横向 测量 标杆 的 间距 为 
10 cm, 纵向 测量 标杆 的 间距 为 20 cm。 通 过 2 台 摄 
像 机 记录 溃 坝 的 全 过 程 。 


(b) 水 档 侧 视 区 


3 试验 设备 及 监测 系统 布设 


Fig.3 Layout of test equipment and monitoring system 


2.2 相似 关系 选择 


土 颗粒 之 间 的 相互 作用 决定 了 土石 坝 尝 口 边 坡 

的 破坏 形式 ( 前 切 破坏 和 张 拉 破坏 ) , 黏 聚 力 c 和 内 

摩擦 角 op 是 影响 颗粒 间 相 互 作 用 的 主要 因素 。 由 于 
砂砾 石料 为 无 笑 性 土 ,选择 模型 材料 时 ,应 尽 可 能 i 

择 c 值 较 小 的 材料 。 内 摩擦 角 w 为 无 量 纲 参数 ,只 

有 确保 模型 内 摩擦 角 与 原型 相等 才能 保证 模型 与 原 

型 省 口 边 坡 的 破坏 形式 相似 

A。 = =1 (10) 

模型 面板 的 选择 除了 满足 长 度 比 尺 Aj 外 ,还 应 

符合 力 的 相似 要 求 。 漫 顶 演 决 过 程 中 ,面板 结构 可 


视 为 悬臂 板 ,主要 受到 3 种 力 的 作用 。 

1) 面板 结构 自身 的 重力 。 重 力 对 面板 结构 的 折 
断 起 着 重要 的 作用 ,为 了 保持 坝 体 与 面板 的 变形 协 
调 性 ,面板 与 砂砾 石料 采用 相同 的 密度 比 尺 A, = 
pa/p, ,由 于 模型 和 原型 中 的 重力 加 速度 相同 ,和 A。= 
1 ,可 以 得 到 重力 比 尺 Ac, 即 

Ac = 和 Ar (11) 

2) 面板 承受 弯 矩 小 于 极限 弯 矩 时 ,产生 弹性 应 
变 的 力 。 根 据 胡 克 定律 ,面板 产生 弹性 变形 的 力 为 弹 
性 应 变 e、 弹 性 模 量 和 截面 面积 4 的 乘积 。 同 时 ， 
为 了 保持 长 度 相似 ,模型 中 的 应 变 必须 等 于 原型 中 的 
应 变 。 因 此 ,产生 应 变 的 弹性 力 比 尺 Ar 被 确定 为 


100 
E 
人 = pa (12) 


基于 力 的 相似 要 求 , 弹 性 力 比 尺 Ar 应 当 等 于 重 

力 比 尺 Ac, 由 式 (11) 和 式 (12) 可 获得 
E 

尿 = Ah， (13) 

3) 面板 折断 时 ,与 面板 极限 抗 拉 强度 相关 的 力 。 
为 了 能 够 合理 地 模拟 面板 江口 演化 规律 ,除了 分 析 
面板 的 受 力 状态 外 ,还 应 当 考 虑 面板 的 强度 特性 , 极 
限 抗 拉 强度 是 面板 破坏 的 主要 特性 ,用 内 聚 力 C 与 
截面 面积 4 的 乘积 来 表示 面板 在 极限 抗 拉 强 度 时 所 
受到 的 力 '*“i。 因 此 ,面板 折断 时 内 聚 力 比 尺 A 为 


Ac = 2 (14) 
志 内 聚 力 比 尺 A6 应 当 等 于 重力 比 尺 和 A, 由 式 
(EB 和 式 (14) 可 获得 
号 = A,AL (15) 


© p 

加 经 过 相似 分 析 ,模型 面板 与 原型 的 极限 抗 拉 强 
庆 绕 值 应 当 与 内 聚 力 比值 相等 ,以 此 对 面板 进行 更 
Wa 的 模拟 。 
2 试验 材料 

原型 坝 最 大 坝 高 70 m, 填 筑 密度 为 2.08 g/cm”， 
内 娠 擦 角 为 40"。 受 水 模 空 间 尺 寸 限 制 , 模 型 坝 坝 高 

>< 
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取 为 0.7m( 长 度 比 尺 A， = 1: 100 ) , 填 筑 料 最 大 粒 
径 为 20 mm。 级 配 曲线 如 图 4 所 示 。 砂 砾石 料 的 密 
度 为 1.86 g/cm (密度 比 尺 和 A,= 0.9; 1 ) ,内 摩擦 角 
为 42"。 为 了 确保 模型 坝 的 均匀 性 和 稳定 性 ,模型 坝 
采用 分 层 铺设 碾 压 ,每 层 铺 料 厚度 为 10 cm, 直 至 填 
筑 到 坝 顶 , 模 型 坝 的 孔隙 率 为 0.45。 

文献 [15 ] 和 文献 [23] 分 别 采用 石膏 面板 、 石 膏 
混合 料 面 板 .仿真 混凝土 面板 进行 模型 试验 ,尽管 这 
些 面板 材料 在 试验 中 是 有 效 的 ,但 均 很 难 满足 密度 
高 和 强度 低 的 特性 。 本 试验 通过 水 泥 、 标 准 砂 重唱 
石粉 (比重 为 4.2 ~4.7 g/cm’, 呈 粉末 状 ,常用 于 加 
重 计 使 用 ) 以 及 水 按 一 定 比 例 制 成 砂浆 面板 材料 。 
根据 力学 性 能 试验 , 表 1 中 的 砂浆 面板 基本 满足 相 
似 要 求 。 


小 于 某 粒 径 百 分 数 /% 


0 
100 10 员 0 0.01 
土 粒 粒 径 /mm 
图 4 面板 砂砾 石 坝 坝 料 级 配 曲 线 
Fig.4 Grading curve of dam fillings for concrete-face 


sand-gravel dam 


表 1 模型 面板 材料 的 配合 比 与 性 能 


© Tab. 1 Proportion and properties of the model face slab material 
和 配合 比 (质量 百分比 ) 性 能 
(水 泥 /% 细 砂 /% 重唱 石粉 /% 水 /% 密度 /(g. em) ” 抗 压强 度 /kPa 弹性 模 量 /GPa ” 抗 拉 强度 /kPa 
10.5 62.6 11.1 15.8 2.13 550 0.28 30.2 


3 数值 模型 设计 


3.1 数值 模型 建立 


为 验证 本 研究 所 建 数学 模型 的 合理 性 ,根据 试 
验 模型 按照 1: 1 建立 数值 计算 模型 。 数 值 模型 坝 高 
为 0.7m, 上 游 坝 坡 1:1.4, 下 游 坝 坡 1:1.5, 坝 顶 宽 
度 6 cm ,面板 宽 10 em, 面板 等 效 厚 度 0.4 em。 另外 ， 
在 坝 顶 中 部 设 有 底 宽 10 cm , 项 宽 14 cm ,深度 为 4 cm 
的 梯形 导 流 槽 。 采 用 结构 化 网 格 对 计算 域 进行 划 
分 ,水 域 采 用 0.04 mx0.04mx0.04m 的 网 格 , 对 


精度 要 求 高 的 面板 和 坝 体 采用 0.02 m x0.02m x 
0. 02 m 网 格 进 行 局 部 加 密 , 网 格 单元 总 数 约 61 万 
个 ,划分 结果 如 图 5 所 示 。 


: 


5 模型 建立 及 网 格 划 分 


Fig.5 Model establishment and mesh division 


3.2 初始 条 件 和 边界 条 件 


根据 模型 试验 的 材料 参数 ,对 砂砾 石料 和 面板 
材料 进行 设 定 。 模 型 的 初始 水 位 高 度 0.57 m, 与 初 
合演 口 底部 高 程 齐 平 。 在 上 游 坡 角 和 坝 顶 处 分 别 布 
设 有 监测 点 和 监测 挡 板 ,以 监测 省 决 过 程 中 库 区 水 
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位 和 溃 口 流量 变化 过 程 线 。 模 型 的 边界 条 件 见 表 
2 ,上 游 为 流量 边界 :人 库 流 量 恒定 为 5 L/s; 下 游 为 
出 流 边界 :流体 从 坝 体 下 游 流出 ,不 发 生 回流 ;模型 
上 部 为 压力 边界 :标准 大 气压 ;模型 空间 的 底部 以 及 
左右 两 侧 为 刚性 墙 边 界 ,相当 于 钢 制 水 槽 的 底部 和 
两 侧 。 


表 2 ， 溃 坝 模型 边界 条 件 


Tab.2 The boundary conditions of the dam-break model 


名 称 Rin i Ys 全 pf 
边界 条 件 刚性 墙 边界 刚性 墙 边界 流量 边界 出 流 边 界 刚性 墙 边界 标准 大 气压 边界 
wl 一 计算 值 
A 一 一 试验 值 
4 计算 结果 分 析 区 

= 0.5 
所 0.4 
水 位 及 流量 结果 对 比分 析 2 
图 6 和 图 7 分 别 是 计算 得 到 的 溃 坝 过 程 中 坝 前 3 
0.1 


双生 溃 口 流量 与 试验 值 的 对 比 ,从 图 中 可 以 看 出 ， 
数值 模拟 较 好 地 反映 了 物理 模型 试验 的 整个 漫 顶 溃 
凑 过 程 ,只 是 在 相位 上 存在 偏差 。 初 始 时 由 于 面板 
站 发 生 破坏 , 随 着 库 水 位 的 抬升 , 溃 口 流量 缓慢 增 
考 叶 流 情 在 下 浊 水 流 的 冲 蚀 与 淘 刷 作用 下 逐渐 展 
宽 和 下 切 ,当面 板 无 法 承受 上 游 水 荷载 和 自重 产生 
的 守 矩 时 ,面板 发 生 第 一 次 折断 。 计 算 模型 在 84 s 
时 爱 生 第 一 次 折断 , 比试 验 值 滞后 了 10 s。 由 于 面 
板 硬 折断 , 溃 口 流量 突然 增 大 ,促进 了 整个 溃 坝 过 程 
的 发 生 ,计算 和 试验 中 均 在 面板 发 生 第 三 次 折断 时 ， 
溃 加 流量 达到 最 大 ,分 别 为 0.068 m*/s 和 0.061 m/s， 
相对 误差 为 9. 8% 。 在 最 大 洪峰 流量 之 后 ,省 口 流 
量 均 在 一 段 时 间 内 保持 较 大 的 同时 进行 衰减 ,最 终 
计算 模型 历时 200 s 演 坝 结束 ,试验 历时 225 s 淡 坝 
结束 ,相对 误差 为 11. 1% 。 分 析 试 验 值 与 计算 值 存 
在 偏差 的 原因 主要 有 两 个 :一 是 , 溃 口 流量 试验 值 是 
通过 坝 体 上 游 坡 脚 处 水 位 计 测 量 出 水 位 ,再 结合 水 
量 平衡 方程 计算 获得 ,而 数值 模拟 是 监测 坝 顶 溃 品 
处 的 流量 ,两 者 溃 口 流量 的 监测 位 置 存在 偏差 ,使 得 
试验 流量 曲线 和 计算 流量 曲线 会 存在 差异 ,但 是 由 
于 两 处 位 置 相距 仅 为 0. 86 m, 对 其 影响 差异 相对 较 
小 ;二 是 ,由 于 砂砾 石料 孔 阶 率 测量 值 偏 大 ,导致 数 
值 模拟 中 坝 体 填 筑 料 比试 验 模型 松散 ,降低 了 坝 体 
填 筑 料 的 抗 冲 蚀 性 能 ,其 数值 计算 的 峰值 流量 偏 大 ， 
省 坝 历 时 偏 短 。 


0 25 50 75 100 125 150 175 200 225 250 
时 间 /s 
到 6 坝 前 水 位 过 程 试验 值 与 计算 值 比较 


Fig.6 Comparison of test and calculated value 


of water level in front of dam 


一 一 计算 值 
0.06 上 一 一 试验 值 


溃 口 流量 /mm ，s 7) 


“0 23 50 73 "100 23 150 175 200 225 250 
时 间 /s 
7 省 口 流量 过 程 试验 值 与 计算 值 比较 


Fig.7 Comparison of test and calculated value 


潭 


of breach discharge hydrographs 
4.2 大 坝 侵 蚀 演 化 过 程 分 析 


面板 砂砾 石 坝 的 整个 漫 顶 演 坝 过 程 分 为 3 个 阶 
段 ,如 图 8 所 示 ( 纵 剖面 )。 第 I 阶 段 :逐渐 冲 蚀 阶段 。 
水 流 流 经 坝 顶 B 点 , 导 流 权 内 出 现 明渠 流 。 在 洪流 作 
用 下 ,细小 颗粒 开始 发 生 运动 ,少量 砂粒 聚集 形成 小 
丘 , 绥 绥 向 下 游 推 移 。 随 着 洪流 量 逐 渐 增 大 ,在 下 游 
坝 顶 4 点 处 出 现 “ 切 头 " 现 象 ,如 图 8(b) 所 示 ,4 点 是 
大 坝 最 先 破 坏 的 地 方 。 当 水 流 达 到 一 定 强度 时 ,在 剪 
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切 力 和 牵引 力 的 作用 下 ,砂粒 从 单 颗 粒 移动 转 为 多 颗 
粒 聚 集 流动 ,下 游 坡 面 逐 渐 冲 出 一 条 冲 沟 ,并 在 坡 脚 
处 形成 鹿 形 堆积 体 。 由 于 面板 参与 挡 水 , 坝 顶 溃 口 处 
形成 瀑布 状 水 流 冲 名 坝 体 ,并 在 落 点 处 冲 出 深 坑 ,如 
图 8(d) 所 示 , 冲 坑 中 的 水 流 为 亲 流 流 态 。 

第 开 阶 段 :剧烈 溃 决 阶段 。 随 着 水 流 的 冲 蚀 , 面 
板 悬 空 长 度 逐 渐 增 大 , 当面 板 无 法 承受 自重 和 上 游 
水 荷载 的 共同 作用 时 ,面板 发 生 第 一 次 折断 (图 8 
(c) ) 。 面 板 折断 使 得 溃 口 流量 突然 增 大 , 冲 坑内 流 
态 的 紊乱 程度 加 剧 ,产生 反 向 湾 涡 ,如 图 8(g) 所 示 。 
游 涡 水 流 一 方面 冲 蚀 下 游 坝 坡 , 另 一 方面 淘 刷 面板 
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下 的 支撑 体 ,加 剧 了 大 坝 的 侵蚀 过 程 。 同 时 ,下 洪水 
流 的 势能 转化 为 动能 ,导致 水 流 以 很 高 的 速度 流 经 
面板 顶部 ,形成 射流 冲 蚀 溃 口 ,如 图 8(h) 所 示 。 随 
着 省 口 的 展 宽 与 下 切 ,面板 支撑 体 陆续 被 侵蚀 ,新 的 
面板 悬空 很 快 形成 ,当面 板 承 受 荷 载 达到 临界 值 时 
再 次 发 生 折断 。 之 后 往复 上 述 过 程 ,直到 水 流 强度 
不 足以 使 面板 发 生 折 断 为 止 。 

第 亚 阶段 :逐渐 稳定 阶段 。 该 阶段 面板 不 再 发 
生 折断 , 库 水 位 逐渐 降低 , 冲 蚀 能 力 减 弱 。 当 省 口 流 
量 趋 于 入 库 流 量 , 床 面 不 再 有 砂粒 的 运 移 时 ,标志 着 


(d) (=90s 


第 I 阶段 ， 逐渐 冲 蚀 阶段 


B 
@ 


第 开 阶 段 ， 剧 烈 溃 决 阶段 


0) {200s 


第 II 阶段 ， 逐渐 稳定 阶段 


图 8 大 坝 溃 决 的 纵 剖 面 侵蚀 演化 过 程 


Fig.8 Longitudinal evolution process of dam breaching 


4.3 溃 口 扩展 过 程 分 析 


图 9 为 试验 与 数值 模拟 过 程 中 部 分 时 段 的 泪 坝 
形态 。 通 过 数值 计算 结果 可 以 得 到 ,第 工 阶段 (0 ~ 
90 s) :由 于 面板 未 发 生 破 坏 , 在 坝 顶 处 形成 跌落 水 
流 ,水 流 在 下 泄 过 程 中 带 走 大 量 的 砂砾 石料 ,逐渐 在 
下 游 坝 坡 上 形成 一 条 沟 字 , 演 口 的 发 展 以 下 切 为 主 ， 
横向 展 宽 为 辅 。 第 开 阶 段 (90 ~ 170 s) :面板 发 生 间 
歇 性 折断 , 溃 口 流量 增 大 使 得 溃 口 不 断 下 切 并 且 溃 
口 边 坡 发 生 失 稳 H 塌 ,该 阶段 江口 的 下 切 和 横向 展 
宽 都 占据 着 主导 作用 。 第 五 阶段 (170 s ~200 s) : 溃 
口 的 下 切 过 程 基 本 结束 , 溃 口 以 横向 展 宽 为 主 。 在 
整个 数值 模拟 中 , 演 口 呈 倒 梯形 ,符合 无 犁 性 土 演 口 


(c) 三 185 s 时 省 坝 形 态 
9 试验 与 数值 模拟 过 程 中 演 坝 形态 
Fig.9 Dam break morphology during test 


and numerical simulation 

图 10 为 坝 顶 溃 口 下 游 1/2 坝 坡 处 省 口 发 展 过 
程 试 验 值 与 计算 值 ,数据 对 比 见 表 3。 计 算 结果 表 
明 ,江口 下 切 深度 与 试验 值 基本 吻合 ,而 江口 横向 扩 
展 与 试验 值 的 相对 误差 均 大 于 15% 。 究 其 原因 , 砂 
砾石 料 虽 为 无 黏 性 土 , 但 在 实际 中 砂砾 石料 仍 是 具 
有 一 定 的 黏 聚 力 , 由 于 在 数值 模型 中 将 砂砾 石料 作 
为 理想 无 黏 性 土 考虑 ,使 得 省 口 横向 扩展 与 试验 中 
溃 口 截面 形状 有 偏差 。 但 从 整体 上 看 ,数值 模拟 和 
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物理 模型 试验 的 结果 基本 一 致 。 


0.7| 一 一 计算 值 (下 游 1/2 坝 坡 ) 
F 游 1/2 坝 坡 ) 


一 一 计算 值 ( 坝 项 ) 
一 一 试验 值 ( 坝 顶 ) 


0.0 
0 25 50 75 100 125 150 175 200 225 250 
时 间 /s 
(a) 溃 口 深度 发 展 过 程 
加 
10F 一 一 计算 值 (下 游 12 坝 坡 ) 
| 一 一 试验 值 (下 游 1/2 坝 坡 ) 
08 2 es 
上 0.7 jf 
其 0.6 
0 
取 04 
一 一 计算 值 ( 坝 顶 ) 
03 一 一 试验 值 ( 坝 顶 ) 


0 25 35075 100" 125 1350 0175 200 225" 230 
时 间 /s 
(b) 江口 宽度 发 展 过 程 
10” 坝 顶 演 口 和 下 游 1/2 坝 坡 处 省 口 发 展 过 程 
Fig. 10 Breach development at the crest and middle 


of the downstream face of the dam 
表 3 江口 发 展 试验 值 与 计算 值 比较 
Tab.3 Comparison of test and calculated 


data of breach development 


参数 计算 值 /m ”试验 值 /m 相对 误差 /( % ) 
坝 顶 演 口 宽度 1.000 0. 820 +22 

坝 顶 淡 口 深度 0.684 0.661 +3.5 

下 游 1/2 坝 坡 溃 口 宽度 。 1.000 0. 851 +17.5 
下 游 1/2 坝 坡 溃 口 深度 。 0.355 0.355 0 


4.4 面板 折断 演化 过 程 


图 11 给 出 了 模型 试验 和 数值 计算 中 6 号 面板 
的 折断 长 度 .折断 时 间 以 及 面板 省 口 的 详细 演化 过 
程 ,计算 数据 与 实测 数据 见 表 4。 由 计算 结果 可 知 , 大 
坝 在 漫 项 84s 后 ,6 号 面板 首次 发 生 折断 ,折断 长 度 为 
12.96 cm ,与 实测 折断 长 度 (14. 30 em) 相对 误差 为 
9.4% ;在 省 坝 170 s 后 ,面板 不 再 发 生 折断 ,6 号 面板 
最 终 垂 向 折断 长 度 为 40. 86 cm ,与 实测 值 (39.35 cm ) 
相对 误差 为 3.8% ,误差 均 控制 在 15% 以 内 。 此 外 , 面 
板 省 口 演化 计算 结果 表明 , 溃 口 中 间 面 板 的 折断 长 度 
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与 折断 时 间 都 呈现 出 先 减 小 后 增 大 的 趋势 ,并 且 两 侧 
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面板 的 极限 折断 长 度 较 中 间 面 板 逐 渐 增 大 。 
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(b) 数值 模拟 中 面板 误 口 演化 过 程 
图 11 试验 与 数值 模拟 中 面板 省 口 演化 过 程 
Fig. 11 Evolution process of face slab breach in the test and numerical simulation 
表 4 面板 折断 长 度 试验 值 与 计算 值 比较 
Tab.4 Comparison of test and calculated breach evolution for each concrete face slab 
面板 ”最 终 垂 向 折断 长 度 ”最 终 垂 向 折断 长 度 相对 误差 / 面板 ”最 终 垂 向 折断 长 度 ”最 终 垂 向 折断 长 度 相对 误差 / 
编号 计算 值 /em 试验 值 /cm % 编号 计算 值 /cm 试验 值 /cm % 
3 24.29 = = 7 39.81 36. 85 +8.0 
4 38.76 17.73 +118.6 8 38.76 31.68 +22.3 
S$ 39.81 32.03 +24.3 9 24.29 21,22 +14.5 
6 40.86 39.35 +3.8 10 0 11.19 —100 


分 析 原 因 是 随 着 省 口 流量 的 增 大 ,省 决 过 程 更 
加 剧烈 ,中 部 面板 承受 的 动 水 压力 增 大 ,导致 面板 的 
折断 长 度 减 小 ; 当 洪 峰 流量 过 后 ,中 部 面板 项 部 的 动 
水 压力 有 所 减 小 ,面板 的 折断 长 度 相 应 增加 ,可 见 水 
荷载 对 中 部 面板 的 折断 起 着 决定 作用 。 而 中 间 面 板 
的 连续 折断 导致 库 水 位 逐渐 降低 ,两 侧面 板 所 承受 
的 水 荷载 大 幅度 降低 , 随 着 溃 口 的 逐渐 扩展 ,两 侧面 
板 悬 空 长 度 增 大 ,在 荷载 的 作用 下 面板 发 生 折断 ,但 
自重 占据 着 主导 作用 ,由 于 面板 的 自重 相对 较 小 ， 
此 越 靠近 演 口 两 侧面 板 的 极限 折断 长 度 越 大 。 面 板 
省 口 演化 计算 绪 果 在 模型 试验 中 得 到 了 验证 ,可 见 
本 研究 采用 弯 矩 平衡 法 能 够 有 效 地 模拟 面板 的 折断 
过 程 。 最 终 整 个 面板 的 折断 演化 呈现 出 锯齿 状 ,并 
且 过 流 断 面 为 倾斜 薄 壁 堰 。 


论 


柄 济 人 研究 成 果 的 基础 上 ,建立 了 一 种 模拟 面板 砂 
硬 硬 坝 漫 顶 溃 坝 的 数学 模型 。 利 用 该 模型 对 面板 砂 
友 答 坝 漫 顶 溃 坝 试验 进行 模拟 计算 ,并 将 计算 结果 


,1 ) 该 模型 着 重 考虑 了 流 涡 水 流 和 砂砾 石料 之 间 
的 病人 蚀 特性 ,通过 RNG hk-s 浇 流 模型 和 流体 体积 法 
(YGF 法 ) 详 细 地 模拟 了 洲 涡 水 流 对 面板 支撑 体 的 
反 员 侵蚀 ,并且 反映 了 溃 决 过 程 中 水 流 特性 、 砂 砾石 
料 刹 混凝土 面板 三 者 之 间 的 相互 作用 ,进一步 的 揭 
表面 板 砂 砾石 坝 的 溃 决 机 理 。 

”2) 基于 弯 和 矩 平衡 法 对 混凝土 面板 的 折断 进行 判 
定 ,计算 了 溃 决 过 程 中 面板 泪 口 的 演化 过 程 。 对 比 
分 析 表 明 ,面板 折断 长 度 和 折断 时 间 的 计算 值 与 试 
验 值 较为 接近 。 此 外 ,省 口中 间 面 板 的 折断 长 度 与 
折断 时 间 都 呈现 出 先 减 小 后 增 大 的 规律 ,并 且 两 侧 
面板 的 极限 折断 长 度 较 中 间 面 板 逐 渐 增 大 ,最 终 整 
个 面板 的 折断 演化 呈现 出 锯齿 状 。 

3 ) 通过 数值 计算 获得 的 溃 口 流量 、 堆 石 体 的 溃 
口 深度 . 溃 坝 历时 以 及 面板 的 折断 长 度 , 与 试验 数据 
之 间 的 相对 误差 均 小 于 15% , 而 堆 石 体 溃 口 的 横向 
扩展 略 有 差异 ,但 从 整体 上 看 ,数值 模拟 和 物理 模型 
试验 的 结果 基本 一 致 。 
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